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Abstract. It is shown that the magnetic current-driven ('kink-type') instability produces flow 
and field patterns with helicity and even with a-effect but only if the magnetic background 
field possesses non-vanishing current helicity B ■ curlB by itself. Fields with positive large-scale 
current helicity lead to negative small-scale kinetic helicity. The resulting a-effect is positive. 
These results are very strict for cylindric setups without z-dependence of the background fields. 
The sign rules also hold for the more complicated cases in spheres where the toroidal fields are 
O ' the result of the action of differential rotation (induced from fossil poloidal fields) at least for 
the case that the global rotation is switched off after the onset of the instability. 
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1. Introduction 



Open questions in stellar physics led to the idea that a dynamo operates in the ra- 
■ diative cores of early- type stars (Spruit 2002). Even the helioseismologic observation of 
rigid rotation of the solar interior shows in this direction. The angular momentum trans- 
port by the large-scale magnetic field pattern (fossil field plus toroidal field induced by 
+3 | differential rotation) does not lead to a solid-body rotation unless the viscosity of the 
plasma exceeds the molecular value by a few orders of magnitude (Riidiger & Kitchatinov 
i— —t\ 1996, Eggenberger et al. 2005). Other examples are given by the evolution of the fast 
rotating early-type stars which can only be understood if i) there is a basic transport of 
angular momentum outwards and ii) the radial mixing of chemicals remains weak (Yoon 
et al. 2006, Brott et al. 2008). Hence, if a (magnetic-induced) instability existed in the 
radiative stellar cores then the corresponding Schmidt number Sc = v j D must be rather 
C***~ ' large. We have shown that the kink-type instability (or Tayler instability, TI) of toroidal 
magnetic fields forms a much-promising candidate for the instability. A Schmidt number 
larger than ten results as the ratio of the effective viscosity and the diffusion coefficient 
; (Riidiger et al. 2009). 

The unstable modes of the TI are basically nonaxisymmetric driven by the energy of the 
electrical current which produces the toroidal field. Interestingly enough, there exists even 
an instability of a toroidal magnetic field which in the fluid is current-free. In this case 
the energy comes from a differential rotation which itself is stable but which is unstable 
under the influence of the (current-free) toroidal field. We have called this instability as 
Azimuthal MagnetoRotational Instability (AMRI) as - like for the standard MRI (with 
axial fields) - the field itself is current-free and does not exert forces. In opposition to 
the standard MRI the AMRI is always nonaxisymmetric and it is, therefore, much more 
interesting for the dynamo theory. For complicated radial profiles of the toroidal field we 
shall always have a mixture of TI and AMRI. Generally, the latter is more important 
for fast rotation ($7 > Ha) and v. v. Here the Alfven frequency J7a for the toroidal field 
is used which derives from the Alfven velocity va — B^/^/JIop as the related frequency. 
Between two cylinders with different radii the toroidal field profile with = AR + B / R 
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(R radius) is free of dissipation. The 'perfect' AMRI appears for A = while the 'perfect' 
TI results for B = 0. 

One can compute the necessary electrical currents to excite both sorts of instabilities 
in a columnar Taylor- Couette experiment with gallium as fluid conductor. As we have 
shown the critical Hartmann numbers for self-excitation of axi- and nonaxisymmetric 
perturbation modes do not depend on the magnetic Prandtl number of the fluid which is 
as small as 10~ 3 for stellar plasma and 10 -5 for liquid sodium (Riidiger & Schultz 2010). 

It is typical for the nonaxisymmetric TI and AMRI that always the two modes with 
m = ±1 are excited for the same critical Hartmann number and also - if supercritical 
- with the same growth rates (Fig. [IJ . Despite their simultaneous existence they can be 
excited as singles with different initial conditions. However, if the initial conditions are 
as neutral as possible with respect to a preferred helicity, in the majority of the cases 
one of the modes dominates after our experiences. 

This behavior may have dramatic consequences with respect to the dynamo theory. 
Both the modes with m = ±1 have opposite helicity with the same total amount. The 
mode with m — — 1 is identical to the mode with m = 1 but in a left-hand system. The 
helicity of m = 1 in the right-hand system equals the helicity of m — — 1 in the left-hand 
system. So it is obvious that in one and the same coordinate system the sum of the 
helicity of m = —1 and m = 1 is zero. As a consequence the instability of a toroidal field 
can only develop helicity if by some reasons one of the modes with m = ±1 dominates the 
other. There are the two possibilities that i) one mode dominates the other by chance (so 
as the matter dominates the antimatter) or ii) the existence of a poloidal field prefers one 
of the modes. We have shown that indeed in stellar radiation zones - if the background 
field has a positive current helicity B ■ curlB - the resulting kinetic helicity (it' • cur hi') 
of the fluctuations is always negative (Gellert et al. 2011). A positive current helicity of 
the background field results if an axial field and an axial electrical current are parallel. 
A negative current helicity of the background field results if an axial field and an axial 
electric current are antiparallel. Hence, the resulting kinetic helicities on the basis of 
current-driven instabilities have, therefore, opposite signs in opposite hemispheres of the 
model. 

2. Cylindric geometry 

We arc interested in the stability of a background field B = (0, B^R), Bo) with Bo = 
const, and the flow u = (0, Rfl(R), 0) in a dissipative conducting fluid rotating between 
two rigid cylinders, v is the kinematic viscosity and rj is the magnetic diffusivity, their 
ratio is the magnetic Prandtl number Pm = v/rj. The stationary background rotation 
law is fl = a + b/R 2 with a and b as constants. fl ln and f2 ou t are the rotation rates of the 
cylinders and B m and B out are the azimuthal magnetic fields there. 

One finds for the current helicity of the background field B ■ curl_B ~ ABq, which may 
be either positive or negative (and of course vanishes for the current-free case A = 0). 

The inner value B m may be normalized with the uniform vertical field, i.e. /3 = B- m / Bq. 
For a profile with B m = B out we have B ■ curl B oc 1//3 for the normalized current helicity 
of the background field. The sign of j3 determines the sign of the current helicity. As usual, 
the toroidal field amplitude is measured by the Hartmann number Ha = B m D / J jilo p v rj 
and the global rotation by the Reynolds number Re = £l ln D 2 /v with D = R ou t — -Rm- 
The Alfven frequency is J1a = B[ n / \JJIop~D. 




The boundary conditions associated with the perturbation equations are no-slip for 
the flows and perfectly conducting for the fields. For all computations it is i? out = 2i?j n . 
For purely toroidal fields the expected net helicities vanish as two nonaxisymmetric 
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Figure 1. The two modes with opposite helicity values of the same amount which can be excited 
with different initial conditions. Re = 0, Ha = 200, /3 = 0, Pm = 1 (from Gellert et al. 2011). 

modes with m = ±1 exist with identical excitation conditions and the same amount of 
helicity but of opposite sign (Fig. [TJ. 
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Figure 2. The normalized negative kinetic helicity and the negative current helicity for the 
nonaxisymmetric perturbations as functions of p. Re = 200, Ha = 100, Q- m = 2f2 ut, Pm = 1. 

We now compute both the helicities for helical fields with various pitch values (3 which 
is inverse to the current helicity of the large-scale field. The astrophysically relevant case 
is > f^A (which is not the classical realization of the Tayler instability) . The instability 
only exists for nonrigid rotation so that we always work with a rotation law oc l/R. 
The radial profile of the toroidal field has been modelled by the above mentioned most 
simple profile with almost uniform toroidal field. Figure [5] gives the results (in units of 
Vt\D). Both the kinetic helicity (u 1 ■ curlu') as well as the current helicity (B' ■ curLB') 
are negative. They are decreasing functions of (3. 

We have also calculated the a-effect via the determination of the electromotive force 
(u 1 x B') by the fluctuations. According to the rule that the azimuthal a-effect is anti- 
correlated with the (kinetic) helicity we expect the azimuthal a-effect as positive for 
f3 > 0. This is indeed the case (see Fig. [31 left). One finds positive and negative signs in 
the container but the positive values dominate so that in the average the azimuthal a- 
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effect is indeed positive. On the other hand, for negative j3 we expect positive small-scale 
kinetic helicity and negative cv-effect. 

The example presented in Fig. prefers to the pitch (3 — 3. For higher values, i.e. smaller 
helicity of the background held, the resulting a-effect becomes smaller and smaller. Figure 
[3] (right) shows that the dynamo number C a = a^D/rj runs as C a ~ C//3 with C of 
order 10~ 2 . It is easy to show that this value is too small by two orders of magnitudes 
to allow the operation of a classical ail-dynamo (Gellert et al. 2011). 




Figure 3. The a-effect for Q > Qa- Left: /3 = 3. Right: The dimensionless dynamo number C, 
of the azimuthal a-effect which decays decays like C//3 with C ~ 0.05. Pm = 1. 



3. Spherical geometry 

We start with the amplification of fossil magnetic fields by shear and the magnetic- 
field back-reaction. The two-dimensional, non-linear simulations thus start with an initial 
differential rotation and a purely poloidal magnetic field in a radiative stellar zone. The 
early phase of the simulation shows a generation and steep amplification of toroidal 
magnetic field through differential rotation. The generated Lorentz force redistributes the 
angular momentum. This is why at the same time of field amplification, the differential 
rotation starts to decrease, and the toroidal-field growth is thus limited. 

Figure [4] shows the maximum magnetic field amplitude in the spherical shell as a func- 
tion of time. The Reynolds number is Re = 20 000 and the magnetic Prandtl number is 
Pm = 1. The dotted and dot-dashed lines in Fig. 2] are the results of a linear stability 
analysis of only the toroidal field. By SI we refer to a velocity perturbation which is 
symmetric with respect to the equator and has m = 1, Al is the corresponding antisym- 
metric perturbation. If at any given time both stability lines are above the solid line, the 
corresponding snapshot is stable against m — 1 perturbations. The stability lines cross 
the solid one at about t = 0.0023 (see Arlt & Riidiger 2011). 

The system is now perturbed in a 3D, nonlinear simulation at a somewhat later time 
(to = 0.003). The perturbation is applied to the magnetic field and has an azimuthal 
wave number of m — 1 (and is symmetric with respect to the equator). The resulting 
flow is thus antisymmetric and can be compared with the Al-mode in the linear stability 
diagram in Fig. |4j 
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Figure 4. The stability limits of each snapshot of the axisymmetric evolution of a sufficiently 
strong initial field. The solid line gives the evolution of the toroidal field. The stability limit is 
given for the nonaxisymmetric perturbation patterns SI and Al (from Arlt & Riidiger 2011). 

The spatial distribution of the resulting helicity is shown in Fig. [5] Note the concentra- 
tion of both the helicities near the inner boundary where the tangent cylinder touches the 
inner sphere. We conclude that a considerable part of the kinetic helicity in the system is 
due to the presence of an inner sphere. This is slightly different for the current helicity. A 
considerable amount of positive current helicity is measured in the bulk of the northern 
hemisphere. The situation is unchanged in runs with perfect conductor boundary condi- 
tions at the inner radius. Interestingly, Reshctnyak (2006) with convection simulations 
also finds negative kinetic helicity along the tangent cylinder. We assume that the nega- 
tive helicity near the tangent cylinder of the northern hemisphere is an inner-boundary 
effect, neither related to convective nor Tayler instability turbulence. 

Figure [5] summarizes the consequences of a numerical experiment. In the unstable 
domain of Fig. [H the instability pattern is calculated without global rotation, i.e. for 
Re = 0. In this case the only existing pseudo-scalar is B ■ cm\B like in the above theory 
in cylindric geometry. The question is whether we can find the same relations between 
large-scale current helicity (left), kinetic helicity (middle) and a-effect (right). Note first 




Figure 5. Left: Kinetic helicity, right: Current helicity. Both after averaging over the 

azimuth. Re = 20 000. 
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that indeed the dominant role of the tangent cylinder for the kinetic helicity and the a- 
effect vanishes so that it is obvious that it is a rotation-induced boundary layer effect. In 
the southern hemisphere the pseudo-scalar B ■ curlB is positive while the kinetic helicity 
is negative (but concentrated to the pole) . The corresponding a-effect given in the right 
panel (and computed also by means of the test-field method) proves to be positive (see 
Arlt & Riidiger 2011). As it must, at the northern hemisphere all signs are opposite. 
Hence, the model fulfills the same sign rules as in the above-discussed cylindric setup. 

4. Summary 

Unstable toroidal fields alone are not able to produce helicity and a-effect. It has been 
shown, however, that helicity and a-effect are produced by unstable magnetic large-scale 
field patterns which themselves possess current helicity B ■ curLB. This is insofar under- 
standable as helicity and a-effect are both pseudo-scalars which can only be nonvanishing 
if in the global setup a nonvanishing large-scale pseudo-scalar like B ■ curlB exists. We 
want to stress, however, that the current helicity of the background field is not the only 
possible pseudo-scalar existing in magnetized stellar radiation zones on which other forms 
of helicity and a-effect may base. 
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Figure 6. Left: B ■ curlB. Middle: kinetic helicity (u' ■ curln'). Right: the azimuthal 
component of the a-effect. Re = 0. 
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